Rupture of a Biomembrane under Dynamic Surface Tension 
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How long a fluid membrane vesicle stressed with a steady ramp of micropipette last before rup- 
ture? Or conversely, how high the surface tension should be to rupture a membrane? To answer 
these challenging questions we have developed a theoretical framework that allows description and 
reproduction of Dynamic Tension Spectroscopy (DTS) observations. The kinetics of the membrane 
rupture under ramps of surface tension is described as a combination of initial pore formation 
followed by Brownian process of the pore radius crossing the time-dependent energy barrier. We 
present the formalism and derive (formal) analytical expression of the survival probability describing 
the fate of the membrane under DTS conditions. Using numerical simulations for the membrane 
prepared in an initial state with a given distribution of times for pore nucleation, we have studied 
the membrane lifetime (or inverse of rupture rate) and distribution of membrane surface tension at 
rupture as a function of membrane characteristics like pore nucleation rate, the energy barrier to 
failure and tension loading rate. It is found that simulations reproduce main features of the experi- 
mental data, particularly, the pore nucleation and pore size diffusion controlled limits of membrane 
rupture dynamics. This approach can also be applied to processes of permeation and pore opening 
in membranes (electroporation, membrane disruption by antimicrobial peptides, vesicle fusion). 
PACS numbers: 05.10.Gg, 87.10.-e, 87.16. A- 



I. INTRODUCTION 

Many aspects of biological life crucially depend on sta- 
bility of cell membranes for which several properties are 
not understood yet. Fluid lipid bilayers are the building 
blocks of biological membranes. Pores in such systems 
play an important role in the diffusion of small molecules 
across biomembrancs Q. As well pore formation is a pos- 
sible mechanism for vesicle fusion Q. In order for any 
vesicle to be useful it must be relatively stable. Yet in 
order to undergo fusion, long-lived holes must occur dur- 
ing the fusion transformation. How membranes actually 
manage to exhibit these two conflicting properties is not 
completely clear, but this can likely be realized only dy- 
namically. Dynamic properties are especially important 
for biological membranes, because their static character- 
istics describe a dead structure whereas life and biological 
functions are associated with molecular motions. Thus, 
the desire to understand the dynamics of biomembrane 
rupture, which is the main aim of this paper, is hardly 
surprising. 

On a microscopic level, pores are formed owing to ther- 
mal motion of lipid molecules and, in principle, various 
types of pores can be distinguished. It usually is as- 
sumed that initially nucleated pores have hydrophobic 
edges; the so-called hydrophobic pores [|| which are spon- 
taneously formed in the lipid matrix. The probability for 
the existence of such hydrophobic pores is determined by 
the free energy of pore as a function of pore radius (see 
Sections [TTI and Hill below). And, when the hydrophobic 
pore exceeds a critical size, a reorientation of the lipids 
takes place converting the pore into an hydrophilic one 
where the head groups form the pore walls. As discussed 



in [4| , these reorientation processes can occur at the very 
early stages of pore formation so that nucleation is the 
crucial step in the rupture process. Note that thermal 
fluctuations can as well lead to transient unstable pores, 
often termed as pre-pores @. In what follows, we will be 
interested on stable pores (i. e., a well defined density of 
pores which can be detected, e.g., by neutron scattering, 
with individual pores forming and resealing revcrsibly) 
and head group reorientation mechanisms will be dis- 
regarded. All details on hydrophobic and hydrophilic 
pores will be lumped into the effective model parameters 
like line tension, (unstressed) membrane surface tension, 
and pore size diffusion coefficient. Electric breakdown 
method provides information on pore size which can be 
drawn from the dependence of membrane conductivity 
on applied voltage [3[ , but dynamical pore characteristics 
can be hardly found by this technique. One of the most 
relevant material parameter controlling pore dynamics is 
membrane surface tension. Surface tension suppressing 
thermal fluctuations and promoting membrane adhesive 
properties, can induce adhesion of a membrane onto a 
substrate or to another membrane, and other tension 
induced morphological transitions including membrane 
rupture. 

A closed vesicle without pores can survive for a very 
long period of time. Pores can form and grow in the fluid- 
lipid membrane in response to thermal fluctuations and 
external influences. Several innovative techniques are 
available for observing transient permeation and open- 
ing of pores. These include mechanical stress, strong 
electric fields (electroporation), optical tweezers, implod- 
ing bubbles, adhesion at a substrate, and puncturing by 
a sharp tip. In all instances, the resulting transient pore 
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is usually unstable and leads to membrane rupture for 
some level of the surface tension. Using the rupturing of 
biomembranes under ramps of surface tension, the chal- 
lenge of the Dynamic Tension Spectroscopy (DTS) is to 
identify and quantify the relevant parameters that govern 
the dynamics of membrane rupture and thereby charac- 
terize the membrane mechanical strength. 




FIG. 1: Image of a 20 mm bilayer vesicle aspirated in a mi- 
cropipette (from Ref. @). 

As a demonstration of the DTS technique, Evans et 
al. @ conducted experiments of rupturing fluid mem- 
brane vesicles with a steady ramp of micropipette suc- 
tion (Fig. [IJ. Rupture tests on different types of giant 
phosphatidychlorine lipid vesicles over loading rate (ten- 
sion/time) from 0.01 and up to 100 mN/m/s produce 
distributions of breakage tensions governed by the ki- 
netic process of membrane failure. One might naively 
expect that lipid membranes to rupture at tensions close 
to hydrocarbon - water surface tension as lipids are held 
together by hydrophobic interactions. However, biomem- 
branes rupture at much lower tension. As pointed out 
by Evans et al., rupture strength of a biomembrane is 
a dynamical property and that the level of strength de- 
pends on the time frame for breakage @- Energy bar- 
riers along the tension driven pathway are determinants 
of membrane strength, and the relative heights of these 
barriers lead to time-dependent changes in strength. To 
describe dynamic tension spectra they observed, Evans et 
al. modeled the membrane breakage as a sequence of two 
successive thermally-activated transitions (dependent of 
loading rates) limited either by specific defect (pre-pore) 
formation or by passage over the cavitation barrier (or 
evolution to an unstable pore) . Accordingly, this descrip- 
tion was formalized into a three-state kinetic model for 
the membrane: the defect-free ground state, the defect 
or metastable state, and the ruptured membrane state 
1- 

Motivated by these experimental and theoretical in- 
vestigations and findings, our objective in this paper is 
to develop a minimal theoretical framework of the DTS 
method to describing the kinetic process of membrane 
breakage. Based on the general framework of Kramers 
reaction rate theory 0, H[, we develop in this paper 
a theoretical framework for DTS to describe the pore 
growth and membrane rupture dynamics as a Markovian 
stochastic process crossing a time-dependent energy bar- 
rier. As mentioned, such a theoretical approach is con- 



ceptually similar to that used by Evans et al., @ (see also, 
[tl [I3l). However, our description is more general than 
that presented in [(| as it characterizes and describes both 
primary nucleation event followed by the continuous dy- 
namics of pore growth and shrinking allowing hence easy 
to follow and adapt for further numerical treatments. 



II. PROBLEM FORMULATION 

In what follows, we treat the membrane as a two- 
dimensional continuum medium and we neglect shape 
fluctuations, i.e., the parameter S = hsT/Aim (Tl| is 
small (e.g., <5 ~ 10~ 3 for lipid bilayers). We will deal 
with thermal fluctuations not related with shape fluctu- 
ations but with the process of barrier crossing for pore 
formation characterized by the parameter e = V^/hsT, 
where Vq is the typical energy cost for pore formation 
at the critical pore radius (see Eq.© below). Thus, our 
investigations will concern the regime, 8 <C lO^ 1 < e. 

Within the framework of the DTS, we describe the 
kinetic of membrane rupture as a succession of two pro- 
cesses: an initial pore nucleation followed by a diffusion 
dynamics of the pore size to membrane rupture. 

A. Pore Formation 

For simplicity, we assume that the net process of pore 
nucleation in a membrane can be described by an acti- 
vated process following a first-order kinetics with a rate 
q, i.e., the distribution of times for the membrane to re- 
main free of pores is given by the exponential distribution 
with the rate, q, which is a function of membrane char- 
acteristics. For the purpose of DTS, we will assume that 
the pore nucleation rate q{cr) is a function of membrane 
surface tension a . 



B. Pore Diffusion 

Once the pore is already formed in the membrane, the 
net energy V(r) of such a membrane of thickness I with 
a circular pore of radius r consists of two opposed terms 
[l2| : the surface tension a, favoring the expansion, and 
the energy cost 7 of forming a pore edge (line tension), 
favoring the closure: 

V(r) = 27T7r — nar 2 . (1) 

Assuming that a > and 7 > 0, and both are constant, 
Eq.([T]) predicts that for r < a, where a = 7/er is the pore 
radius for the maximum energy V(r), the radial force as- 
sociated with a change in radius tends to reseal the pore, 
and the membrane remains stable against pore growth. 
On the other hand, a pore with a radius larger than the 
critical value a will grow without bound and, ultimately, 
will rupture the membrane. In DTS experiments @, the 
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membrane is stressed such that (provided that 7 remains 
constant) the surface tension grows linearly with time as, 
a = <tq + Ft, where <tq is the unstressed membrane ten- 
sion and F is the loading rate constant. In this case, the 
critical radius a(t) becomes a decreasing function of time 
and any pore initially with radius r < a(0) will ultimately 
lead to membrane rupture at time such that r > a(t) as 
a result of the decreasing of both the critical pore ra- 
dius and associated barrier energy. Now, incorporating 
thermal fluctuations in this picture, one can view the 
rupture of the membrane as a Brownian process crossing 
the time-dependent energy barrier V[a(i)]. 

To setup the equations of pore size dynamics, we con- 
sider a membrane with a pore of radius r under mechan- 
ical stress that changes its surface tension. In this de- 
scription model of DTS, the surface tension a in V(r) is 
a linear function of time as defined above. Thus, neglect- 
ing incrtial effects, the dynamics of r is governed by the 
Langevin equation with time-dependent potential, 



k V(r, t) = 2tt7 r - 7T (cr + Ft) r 2 



(2) 



where £ = 47r?/ m I is the friction coefficient to radial circu- 
lar fluctuations with rj m the internal 2d membrane viscos- 
ity, and f(t) is a Gaussian random force of zero mean with 
correlation function given by the fluctuation-dissipation 
relation, (f{t)f{t')) = 2(k B TS{t - t'), with k B T being 
the thermal energy. 



TABLE I: Parameters and dimensionless variables. The index 
"0" denotes quantities for unstressed membrane. 



Symbols 


Definition 
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line tension (energy /length) 


on 




unstressed surface tension (energy /surf ace) 


F 




tension loading rate (energy /surf ace/time) 


D = 


-- k B T/Cl 


pore diffusion coefficient (length 2 /time) 


ro - 


= l/aa 


critical pore radius (length) 


T — 


rl/D 


diffusing time scale of the critical pore (time) 


F 


= <?o/t 


critical tension loading rate (energy /surface/time) 


'/o 




reduced unstressed pore nuclcation rate 


x — 


r/r 


reduced pore radius 


y = 


cr/a 


reduced membrane surface tension 


e = 




reduced energy barrier for unstressed membrane 


V — 


F/F 


, 1 , . , I v < 1 : diffusing limit 
reduced tension loading ratc:< 




[ v > 1 : drift limit 



and we have defined the potential (energy landscape for 
DTS illustrated in FigdJ), 



U(x\y) = - 



= - \2x 



yx 



V(ro) 
k B T <7ok B T 



7T7 



(5) 



The potential U(x\y) is maximum at x$ = l/y corre- 
sponding to the energy barrier = U(x*\y) = s/2y. 
Both the position and height of the energy bar- 
rier decrease as y gets larger as a result of the membrane 
stress. 



C. Dimensionless equations 

To work with dimensionless variables, we define in Ta- 
ble U scales of length, surface tension and time by r , <r 
and t, respectively, and we consider the following trans- 
formations: x = r/Vo, y = ct/cto, and t —> tjr (with 
x e [0,1] and y G [l,oo[) . This operation leads us to 
define the control parameter, 



Ft _ diffusing time scale 
do surface tension time scale 



(3) 



This parameter allows us to distinguish two regimes in 
the dynamics of the membrane rupture: the diffusion 
controlled regime when v <C 1 and the drift regime for 
v ^> 1 limit. With these transformations Eq.@ can be 
rewritten as, 



dx 
~dt 

dy 
dt 



dU(x,y) 
dx 



X(t) 



(4) 




tension, y 



pore radius, x 



FIG. 2: Energy landscape (Eq.© with e = 2) for the Dy- 
namic Tension Spectroscopy. Flat area x > 2/y corresponds 
to U < 0. 



D. DTS observables 



where X{t) is a Gaussian random force of zero mean with 
correlation function given by, (X{t)X(t')) = 2 5(t — t'), 



As the barrier crossing to both pore nuclcation and 
membrane rupture are stochastic processes, both the 



4 



membrane lifetime and the membrane tension at rupture 
are distributed. Our goal is to calculate the two quanti- 
ties that characterize the kinetics of membrane rupture 
in DTS framework: the rate of membrane rupture and 
the distribution of tension at membrane rupture. 

For a membrane with a pore in the absence of mechan- 
ical stress v = 0. the distribution of tension at membrane 
rupture is delta function, Q(y) = S(y — 1), and the rate 
of membrane rupture can be obtained analytically using 
the first passage time approach [IH, [l4[ , 

c 1 2 

p eq (z\l)dz 



1 



fc(e|0) 



dx 



Pcq(x|l) 

V^7(2i) 



crfi 



G/2 



dx e~ 



e{x-l) 2 /{2) 



crfi 



— erfi 



{l-x)^7j2\} 



(6) 



in which we have assumed that the membrane system was 
initially prepared with the distribution p eq (x|l), where 



p cq (x\y) 

and Z(y) = 



2ey 



; Z(y) -- 

1/2 

crfi 



dx , (7) 



yje/2y where erfi[z] 



crf[iz]/i and erf[- ■ • ] is the error function. 

On the other hand, for a membrane initially free of pore 
and for v > 0, analytical expressions are not straightfor- 
ward but the rate of membrane rupture and the distribu- 
tion of tension at membrane rupture can be determined 
as follows. Let S(t) be the survival probability that de- 
scribes the fate of the membrane from the beginning of 
the experiment. The distribution of membrane lifetime 
or rupture time is given by —dS/dt, and the membrane 
rupture rate (equals to the inverse of the membrane life- 
time) is obtained as, 



S{t) dt . (8) 



k(e\v) 



poo 

/ 1 


' dS{t)~ 


dt= [ 


10 


dt 


Jo 



Likewise, the distribution Q(y) of tensions y at which 
the membrane rupture is related to the distribution of 
rupture time and, as y = 1 + vt, we have: 



Q(y) 



dS_ 

' dt 



(9) 



y=l+vt 

The DTS spectrum (here, mean of rupture tensions), 

/oo 
yQ(y)dy, is related to the rupture rate 

k(e\v) by, (y(e\v)) = 1 + v/k(e\v). In the case where 
S(t) satisfies a first-order rate equation with the time- 
dependent rate F(<), i. e., S(t) = exp j-^ r(t')eft'j, 
the distribution Q(y) can written as, 

r(y-l)/v } 

T(z)dz } . (10) 



r(2^i) 
Q(y) = v 2 ; exp. 



It follows from this that S(t) is the key function to deriva- 
tion of expressions of quantities of interest. 



III. ANALYTICAL THEORY 

Equivalently to the stochastic equation in Eq.Q, the 
joint probability density, P(x,y,t), of finding the mem- 
brane with surface tension y and pore of radius x (i.e., 
the phase space point (x, y)) at time t is described by the 
two-dimensional Fokkcr-Planck equation: 



dP(x,y,t) 
dt 



dP{x,y,t) dJ(x,y,t) 



dx 



(11) 



where the first term in the right hand side describes the 
ballistic drift of the surface tension caused by the applied 
loading rate, and the second term is the diffusive flux 
describing the diffusion of the pore radius in the potential 
U(x\y) for given y, 

J(x, y,t) = - e~ u ^ ^- e u ^ P(x, y, i) . (12) 
ox 

Eg. pip reduces to the Smoluchowski equation in the 
v = limit 15(. To study the rupture of membrane as 
an escape of the pore radius undergoing a Brownian dy- 
namics within the interval [0, 1 [, we require that P(x,y,t) 
satisfies the reflecting boundary condition at x — and 
the absorbing boundary condition at x = 1/y, i.e.: 



J(x,y,t)=0 at x = , 
P(x,y,t)—0 at x = l/y. 



(13) 



Formal, yet numerically computable, solution of Eq. (|ll[) 
with the initial condition P(x,y,t = to|^0;J/o) = S(x — 
x o) $(y ~ Ho) an d boundary conditions in Eq. (fT3")) is given 
by the Green's function, 



P{x,y,t\x ,y ,t ) 



Peq(x\y) 



-,1/2 



Peq{x \y ) 



S[y-yo- v(t - 1 )] 



/JV'n(^obo) ipn(x\y) exp 



\ n (z) dz 



(14) 



The i\) n (x\y) and X n {y) are respectively the normalized 

( f 1/V \ 
eigenfunctions / dxtp n (x)il ! n'(x) = S n , n ' and asso- 



ciated eigenvalues of the eigenvalue problem, 
" 2 e 2 (l — yx) 2 



Hi> - 

dx 2 



vex 



ey 
2 



i> = -AV> (15) 



satisfying the reflecting and absorbing boundary condi- 
tions at x — and x — 1/y, respectively, 



-U(x) 



d r 



JJ{x)/2 



dx 

. tp(x = 1/y) 



ip(x) 







(16) 
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Let z = [ve + e 2 y 2 ] 
Z\, where, 



1/4 



ey 



v + ey 2 



such that zq < z < 



ey [ve + e 2 y 2 ] 



1/4 



Zl 



v [ve + e 2 y 2 ] 



1/4 



(v + ey 2 ) y{v + ey 2 ) 

The eigenvalue problem, Hip = —Xip, becomes, 

..2- 



(17) 



d 2 ^ 
dx 2 

E 



E 



z 

T 



i 



[ve + e 2 y' 



,1/2 



i' = o 

A - 



ey_ 
2 



ve 



A(v + ey 2 ) 



(18) 



The general solution to Eq. (|T8]) which satisfies the ab- 
sorbing boundary condition in Eq. (|16[) is given by, 



tl>(z) = A [D v {- Zl )D v {z) - D v { Zl )D v {-z)] 
v = E-- 

where D v {z) is the Weber's parabolic cylinder function 
[l6j. The constant A is obtained from the normaliza- 
tion and the eigenvalue spectrum by solving the follow- 
ing eigenvalue equation obtained be using the reflecting 
boundary in Eq. pi 



-zi. 



D v ( Zl ) 
= 0. 



dD v (z ) £ 

dzo 2 [ve + e 
dD v (-z ) e 



V] V4 



dzn 



2[ve + e 2 y 2 } 1/4 



DM 



D v (-z ) 



(20) 



Now, assuming that the system is initially prepared 
with the distribution g{x, y, t) describing the pore for- 
mation, the survival probability that describes the fate 
of the membrane with a pore is given by, 



S(t) 



r°° rl/yo rl/y 
dyo / dy dx / dx 

Ji Jo Jo 

P(x, y, t\x , y , t ) g(x Q ,y , t ) , 



die 



(21) 



where P(x,y,t) is given in Eq.(|14p and the preparation 
distribution, 



g{x,y,t) = p eq (x\y)5(y- 1 -vt) 



q(t) exp 



q{t')dt' 



(22) 



where the term between squared brackets stands for the 
distribution of times for pore nucleation. 

Equation (f2~Tj) provides an exact expression of S(t) in 
terms of infinite series from which the rupture rate k(e\v) 
(or the DTS spectrum) and the distribution Q(y) of rup- 
ture tension can subsequently derived by using Eqs.© 
and dH]), respectively, and approximate expression for T(t) 
as well. 



Interestingly, these derivations can also be used to es- 
tablish the correspondence with the three-state model in 
Ref. @ and, therefore, provide exact expressions as, 



S °( t ) = Sr d V fo /V ° dx ° Io dt ° 

P(0, y, t\x , y , h) g(x , y , t ) 
S.(t) = S(t)-S (t) 

Shole(t) = l-S(t) 



(23) 



where (t) is the defect-free ground state, S 1 * (t) the de- 
fect or metastable state, and SW,i e (£) the ruptured mem- 
brane state as defined in Ref. [6[. 

Unfortunately, derivation of analytical expressions 
(which require solving Eq. (|20[0 as outlined above may 
be tedious and obtained results turn out not easy to use 
in practice. These calculations were done mainly for the 
purpose of presenting the derivation formalism of exact 
expressions. Such exact solutions may turn out useful 
for checking simulation results just like those presented 
in the next section. Our main interest in this paper is 
to understand, write down equations describing DTS ex- 
periments and develop related simulations that could be 
compared with experimental data. To this end, we switch 
to the simulations of the kinetics of the membrane rup- 
ture as described by stochastic and dynamical equations 
outlined above. 



IV. SIMULATION ALGORITHM 

The simulations of the kinetics of the membrane rup- 
ture were performed using the discretized version of 
Eq.@ to have the algorithm, 



x n+ i = {eAy n + l)x n - eA + X n 
Vn+l = y n + vA 



(24) 



where A is the time step and the Gaussian random 
noise X n is defined by the moments, (X n ) = and 
(X n X n i) = 2A6 nn /. For each trajectory for a mem- 
brane free of pore at t = 0, with fixed barrier height 
e and loading rate v, the system is prepared according to 
the distribution g(xo,yo,to) given in Eq. (|22[) : the initial 
pore is created in the membrane at time to = (jjo — 1)A' 
where the membrane surface tension yo for pore creation 
is given by the exponential distribution, 



f(y) = l qM ) 



dy' 



(25) 



and the pore radius xq is generated from the distribution 
Peq(xo\yo) in Eq.([7]). From this, the next pore radii x„ 
and surface tensions y n are generated according to the al- 
gorithm in Eq. (|24p . To simulate the rupture of the mem- 
brane, each trajectory starting at xq (0 < xq < 1/yo) 
at time t = t is terminated at time ^ = nA when the 
condition x n > l/y n is satisfied for the first time (the 
boundary at x — is reflecting). The rupture surface 
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tension j/j = y n , the first passage time U and the sur- 
vival probability Si(t) (defined as Si(t) = 1 for all t < ti 
and Si(t) = otherwise) for this given trajectory are 
recorded. The distribution Q(y) of rupture tensions is 
obtained by binning the j/j's over a large number N of 
trajectories. Likewise, the definitive survival probability, 
S(t), and the rupture rate constant, k{e\v), (i.e., the in- 
verse of the membrane mean lifetime) are then obtained 
by averaging these quantities over a large number of tra- 
jectories: 

N N 

S[t) = -Y,Sm and = J>. (26) 

i=l y ' ' i=l 

For all simulations reported in this paper we used A = 
10~ 5 and a total of N = 10 5 trajectories were used to 
perform the averages. 



V. RESULTS 



B. Finite Pore Nucleation Rate q limit 

To investigate the effect of pore nucleation rate on the 
membrane rupture, we consider two cases of increasing 
complexity. 

1. a — limit 

In this case, fc n ,eff = = Qo- Figure [5] shows the 
variation of rupture rate k(e\v) as a function of the pore 
nucleation rate qo- As expected for a successive process, 
k(e\v) linearly grows with go ui the nucleation controlled 
limit when go *S kd and saturates to kd in the diffusion 
controlled limit for q ;» kd- Figure [6] illustrates that 
system parameters can be tuned to follow the transition 
between the nucleation and diffusion controlled limits. 
Accordingly, Fig. [7] displays the profiles of the distribu- 
tion Q(y) corresponding to the nucleation controlled limit 
and toward the diffusion controlled limit. 



In what follows, simulations were carried out with 
the tension-dependent pore nucleation rate given by, 
q(y) = qoe"^ -1 ' , where qo is the pore creation rate in 
the unstressed membrane and a is a constant depending 
on membrane characteristics and temperature (e.g., s). 
As the membrane tension increases with time with load 
v the overall rate of pore formation will be given by, 



Qo 



dx exp 



1) 



Qo 



(27) 



The DTS outputs are the distribution Q(y) of tensions 
at rupture and the DTS spectrum defined by the plot 
of the mode of Q(y) as a function of ln(v) [§[. We have 
computed the membrane survival probability (results not 
reported) distribution Q(y) and rupture rate k{e\v). As 
a successive process, the membrane rupture rate can be 
written as, k(e\v) = fc d fc, liCff / (k d + fc n ,eff) where fc njCff 
(different from k n ) is the effective rate of pore forma- 
tion and kd is the rupture rate for a membrane initially 
with a pore in it. In what follows, we will investigate the 
effect of q, v and e and on Q(y) and k{e\v). 



Diffusion Controlled Limit: 



oo limit 



This limit corresponds to the situation where mem- 
brane under tension stress has already a pore in it and, 
therefore, k(e\v) = kd- Figures [3l [4] and [5] display the 
distribution Q(y) and the DTS spectrum as a function 
of the energy barrier e and loading tension rate v. As 
can be seen, the distribution Q(y) of tensions at rupture 
broaden from the delta function at y = 1 to a wider dis- 
tribution when both v and e increase. Accordingly, the 
mean tensions (y(e\v)) for membrane rupture increase 
with both loading rate v and barrier height e. 



2. Case of a — 1 

When a ^ 0, the dynamic nature of membrane rupture 
leads to fc n ,eff ^ k n . Both k(e\v) and Q(y) still exhibit 
similar behaviors observed in the case of a = but with 
nontrivial dependence on the loading rate v and barrier 
height e. As kd(s\v) is known from the limit q —> oo, 
the behavior of k nyB s(e\v) can be learned from Fig. [5] 
Clearly, fc njC ff > fc n , and the effective pore nucleation 
rate has extra v and e dependencies that are not taken 
into account in k n in the absence of pore dynamics. The 
departure of k n ^ c g from k n increases with both v and 
e indicating that opening pore is more likely for high 
barrier membrane with high loading rates. The kind of 
distributions Q(y) that can be observed in this limit are 
displayed in Fig. [SJ 



VI. CONCLUDING REMARKS 

Rupturing fluid membranes or vesicles with a steady 
ramp of micropipette suction produces a distribution of 
breakage tensions governed by the kinetic process. Ex- 
perimental evidences have demonstrated that the mem- 
brane rupture is a dynamical property whose the strength 
depends on the time scale for breakage. From the theo- 
retical point of view, we have developed a minimal model 
for the Dynamic Tension Spectroscopy that describes the 
pore nucleation as a first-order activated process, the dy- 
namics of pore growth as a two-dimensional (in pore ra- 
dius and tension spaces) Markovian stochastic process, 
and the rupturing of membrane is modeled by an escape 
process over the time-dependent critical barrier of the 
energy landscape. We have provided an exact analytical 
solution of this problem and established the correspon- 
dence between this description and the three-state model 



7 



in Rcf. [6|. As numerical results, we have simulated the 
rupture rate and the distribution of rupture tension as 
a function of the pore nucleation rate, the critical bar- 
rier height, e, of the unstressed membrane and the re- 
duced tension loading rate, v. Our simulated histograms 
reproduce already several features observed in DTS ex- 
periments in [(| and highlight the variety of profiles and 
richness of the problem. Indeed, the distribution of rup- 
ture tensions show different profiles between and in the 
two nucleation and diffusion controlled limits as a func- 
tion of e and v. 

As presented above, the kinetic of membrane rupture 
as probed in DTS experiments is very similar to non- 
equilibrium problems studied in single-molecule pulling 
experiments using atomic-force microscopes [17 1 llal . To 
cite a few, there are several theoretical works [19l - |23l ] that 
have been developed, extended and refined to describe 
the thermally activated rupture events within the general 



framework of Kramers reaction rate theory [7j, |8[ . 

Needless to recall that the theoretical model outlined 
above does not yet take into account all aspects of the 
membrane rupture observable in the experiments since 
we purposely neglected some of features like, e.g., the 
coarse-grained structure of the membrane. However, the 
developed formulation can be embellished in several di- 
rections to include the mentioned above and some other 
ingredients, like, non-Markovian dynamics of the pore ra- 
dius dynamics driven by the membrane matrix in which 
the pore is embedded, and eventually, the time varia- 
tion of barrier energy e due to change of the line tension. 
Such a generalized reaction rate approach can be also ap- 
plied to membrane disruption by antimicrobial peptides. 
Indeed a peptide binding causes a local membrane area 
expansion and therefore it is equivalent to local tension 
(see more about the problem in 0, [l^ ) . Such a work is 
underway. 
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FIG. 3: Distribution, Q(y), of rupture tensions, y, for vari- 
ous values of the energy barrier, e, and loading rate, v. For 
each row of figures, the leftmost and second leftmost figures 
have the same scales in both x and y axes, and likewise for 
rightmost and second rightmost figures. Note that y - scales 
of rightmost figures are an order of magnitude smaller than 
for the leftmost ones. 
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FIG. 4: DTS spectrum (y(e\v)) as a function of the load- 
ing rates v (in logarithm scale) for various barrier heights e 
(quoted numbers). Filled symbols correspond to simulations 
and dashed lines are guide eyes. 




FIG. 5: DTS spectrum (y(e\v)) as a function of the barrier 
height e for various loading rates v (quoted numbers). Filled 
symbols correspond to simulations and dashed lines are guide 
eyes. 
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FIG. 6: Rupture rate k(e\v) as a function of the tension inde- 
pendent pore nucleation rate go (for a — 0) for various barrier 
heights e (quoted numbers) and loading rates v = 0.1 (circles 
and squares) and v = 10 (diamonds and triangles). Filled 
symbols correspond to simulations and dashed lines are guide 
eyes. 
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FIG. 7: Distribution, Q(y), of rupture tensions, y, for various 
values of the tension independent pore nucleation rate, qo (for 
a — 0), and loading rate, v. For each row of figures, panels 
with qo = 1 and go = 10 have the same scales in both x and 
y axes. 
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FIG. 8: Reduced effective nucleation rate k n ^a{e\v) / k n {v) as 
a function of the loading rates v in log-linear scale for various 
barrier heights e (quoted numbers), with go = 0.1 and a = 1. 
Filled symbols correspond to simulations and dashed lines are 
guide eyes. Inset: Reduced nucleation rate k n /qo versus the 
variable av/qo in log- log scale. 




FIG. 9: Distribution, Q(y), of rupture tensions, y, for various 
values of loading rate, v, and energy barrier, e, with qo — 0.1 
and a = 1. For each row of figures, panels with v = 1, v = 10 
and w = 100 have the same scales in both x and y axes. 



